Multiscaling in Hall-Magnetohydrodynamic Turbulence: Insights from a Shell Model 
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We show that a shell- model version of the three-dimensional Hall-magnetohydrodynamic (3D Hall- 
MHD) equations provides a natural theoretical model for investigating the multiscaling behaviors 
of velocity and magnetic structure functions. We carry out extensive numerical studies of this shell 
model, obtain the scaling exponents for its structure functions, in both the low-fc and high-fc power- 
law ranges of 3D Hall-MHD, and find that the extended-self-similarity (ESS) procedure is helpful 
in extracting the multiscaling nature of structure functions in the high-fc regime, which otherwise 
appears to display simple scaling. Our results shed light on intriguing solar-wind measurements. 
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Turbulent plasmas abound in accretion disks, galaxies, 
stars, the solar wind, and laboratory experiments [2 12]; 
thus, the characterization of the statistical properties [1, 
[31 H] of turbulence in such plasmas is a problem of cen- 
tral importance in astrophysics, plasma physics, fluid dy- 
namics, and nonequilibrium statistical mechanics. Such a 
characterization begins with the energy spectra: e.g., in 
homogeneous and isotropic fluid turbulence the energy 
spectrum E(k), which gives the distribution of energy 
over different wave numbers k, assumes the scaling form 
E(k) ~ k~ a if the Reynolds numbers Re is large and 
k is in the inertial range L^ 1 <C k <C k d , where L is 
the energy-injection length scale and kd = 2ir/rjd, with 
i]d the length scale at which viscous dissipation becomes 
significant; the phenomcnological theory of Kolmogorov 
(K41) yields H] a — 5/3. Turbulent plasmas show 
similar scaling forms for the kinetic and magnetic-energy 
spectra E u (k) and E b (k), if the turbulence is statisti- 
cally homogeneous and isotropic, and both Re and the 
magnetic Reynolds numbers ReM are large; their ratio 
PrM = ReM /Re, the magnetic Prandtl number, governs 
the relative sizes of the fluid and magnetic dissipation 
length scales rf^ and rj d ; the inertial-range scaling proper- 
ties of E u {k) and E b (k) have been studied theoretically 
and numerically by using the equations of magnetohy- 
drodynamics (MHD) [2I31H]- Energy-spectra measure- 
ments in the solar wind [7] have shown, however, that 
E b (k) displays two power-law ranges. Several authors j8j- 
[P2] have suggested that, to obtain these two power-law 
regimes, we must augment the MHD equations with a 
Hall-effect term, which leads to a scale separation at the 
ion-inertial length dj or, equivalently, at the wave number 
ki = l^jdi- For k < kj, E b (k) oc k~ a ' , it has been ob- 
served that a M ~ 5/3. For k d > k > k T) E b (k) cx k- a "' 2 , 
where a b ' 2 is either ~ 7/3 or ~ 11/3. The value of a b ' 2 
depends on whether the magnetic energy dominates over 
the fluid kinetic energy, which occurs in the electron- 



MHD (EMHD) [13] limit, or the converse, i.e., the ion- 
MHD (IMHD) limit. These limits follow from the 3D 
Hall-MHD equations: EMHD is obtained if the induc- 
tion term is sub-dominant to the Hall term; in the IMHD 
case these two terms are comparable to each other. In the 
EMHD limit, we obtain a single, characteristic scale and 
K41 phenomenology yields a b ' 2 = 7/3; in the IMHD case 
a comparison of the transfer time, from the Hall-term, 
and a second time, from the induction part, followed by 
simple dimensional analysis yields a b:2 = 11/3 [TUI IT2] . 

Direct numerical simulations (DNSs) [9hTT| have just 
begun to resolve these two scaling ranges; but their 
spatial resolution is much more limited than has been 
achieved in DNS studies of MHD turbulence 0H]. Thus, 
they have not been used to study the scaling or multi- 
scaling properties of order p fluid and magnetic struc- 
ture functions (defined below). However, measurements 
of such equal-time magnetic structure functions in solar- 
wind measurements [7] show that, although there is sig- 
nificant multiscaling in the low-fc {k < kj), power-law 
range of E b (k), the scaling exponents in the second, high- 
k (kd > k > kj) power-law range increase linearly with 
the order p. Thus, it behooves us to develop a theo- 
retical understanding of these important and intriguing 
observations and to test them. 

We show that a shell-model version of the 3D Hall- 
MHD equations [TTJ [T2], which is a generalization of 
MHD shell models O [TS] , provides a natural theoreti- 
cal model for investigating such multiscaling behaviors in 
structure functions in 3D Hall-MHD turbulence. Given 
the large range of scales that we can cover in this shell 
model [TBI, its magnetic spectrum E b (k) reveals two, dis- 
tinct, power-law ranges. We carry out the most compre- 
hensive numerical study of this 3D Hall-MHD shell model 
attempted so far; and thereby we characterize and quan- 
tify, for the first time, the properties of the order-p mag- 
netic and velocity structure functions in this model via 
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10" 8 


5 x 10" 13 
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0.1 


1.2 ±0.3 


1.8 ±0.2 


1.6± 0.4 
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5 x 10~ 12 
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0.1 


4.3 ±0.8 


4.7 ±0.4 


1.1 ±0.2 
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1 x 10" 13 


1 x 10~ 14 


0.01 


2.8 ±0.6 


0.74 ±0.09 


0.27 ±0.06 
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1 x KT 14 


0.01 


1.0 ±0.2 


0.31 ± 0.04 


0.31 ±0.07 
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a b,2 


Rl 


16.3 ±0.6 


6.64 ±0.09 


11 ±2 


1.69 ±0.04 


2.2 ±0.3 


1.66 ±0.03 


1.75 ±0.03 


3.39 ±0.08 


R2 


16 ±1 


6.7 ±0.1 


4.9 ±0.8 


3.40 ±0.08 


1.7 ±0.2 


1.65 ±0.03 


1.69 ±0.03 


3.45 ±0.06 


R3 


16 ± 1 


6.6 ±0.1 


26 ±4 


0.90 ±0.03 


0.62 ±0.07 


1.69 ±0.03 


1.74 ± 0.02 


3.34 ±0.09 


R4 


13.3 ±0.7 


6.8 ±0.1 


64 ±9 


0.46 ±0.02 


0.64 ±0.07 


1.70 ±0.02 


1.74 ± 0.02 


3.28 ±0.08 



TABLE I: The values of the different parameters (see text) used in our runs Rl — R4 and the spectral exponents a", a 6 ' 1 , and 
a b ' 2 obtained. 



their scaling exponents Cp (fluid), (p 1 (magnetic, k < ki 
regime), and Q b ' 2 (magnetic, kd > k > fcj regime). We 
find that all three sets of exponents show clear signatures 
of multiscaling. In particular, we find the remarkable re- 
sult that magnetic structure functions display multiscal- 
ing for both the low-fc and the high-/c power-law ranges. A 
second significant and surprising finding is that, although 
the exponents £p 2 ^ ( b ' x , ( b ' 2 ^ Q, the exponent ratios 

The 3D Hall-MHD equations for the velocity u and 
magnetic b fields are 

-t^±(u-V)u = -Vp+jxb + i/V 2 u; 

°— = Vx [(u-d;j)xb]+,,V 2 b; (1) 

here v and 77 are the kinematic viscosity and magnetic 
diffusivity, respectively, di is the ion-inertial length, the 
scale at which the Hall effect becomes important, the cur- 
rent density vector j = V x b, the pressure is p, V ■ b = O, 
and, at low Mach numbers, the flow is incompressible, 
i.e., V ■ u = O. We define the dissipation length scales 
7$ = (i^/e") 1 ^ and ^ = (r/ 3 /£ b ) 1/4 , where e u and e b 
are the kinetic and magnetic-energy dissipation rates, re- 
spectively; we restrict ourselves to decaying turbulence, 
so we do not include forcing terms. The Hall term, which 
is a singular perturbation of the MHD equations [TU] , has 
a significant effect if di ^> n^, rj h d . The shell- model ver- 
sions of Eq.Q are PUCE]: 

— j = -vk 2 n u n - v 2 k^u n + t^]*, 

^ = -nklbn-^kX + (2) 

where u n and b n are, respectively, the complex ve- 
locity and magnetic field in the shell n, * denotes 
complex conjugation, 1 < n < N, where N is the total 
number of shells, $JJ = A n {u n+1 u n+2 - b n+1 b n+2 ) + 



B n (u n -lU n+ i — b n -lb n+ i) ± C n (u n - 2 Un-l — & n _ 2 &ri-l) 

and $ b n = D n (u n+ ib n+2 - b n+ iu n+2 ) + E n (u„-ib n+ i - 
b n -iUn+i)-Fn(u n -2bn-i~b n _ 2 u n _i)~d I [G n b n+ ib n+2 + 
H n b n -iK+i ± 7„6„_2&n-i], with A n = k„, 
Bn 2 kn—ii C n 2 k n — 2j D n Qk rlJ E n — ^k n — 1, 

F — lb n H — -if— IV'P H — — iC — 1 V^P 

I n = (-l) n ~ 2 k n _ 2 , k n = 2 n k , and k = 1/16; the values 
of the coefficients A n — I n are determined by enforcing 
the shell-model analogs of the Hall-MHD conservation 
laws, in the inviscid, unforced limit; the conserved quan- 
tities are the total energy E = E„(|w„| 2 ± |&„| 2 )/2, the 
magnetic hclicity Hm = S„(- l) n \b n \ 2 /2k n , and the ion 
helicity Hi = E n ((&„< + b* n u n ) + di(-l) n k n \u n \ 2 /2); 
the hyperviscosity v 2 and the magnetic hyperdif- 
fusivity r\ 2 have to be included for numerical sta- 
bility [TQl [T3]. We use the boundary conditions 
An-i = An — B\ = Bn = C\ = C 2 = 0, 
Dn-i — Dn = Ei — En = i*i = F 2 = 0, 
Gjv_i = Gn = Hi = H N = Ii = Ii = 0, and the 
following initial values for u n = u^°' ) k n 1 ^ 3 e^ k " +l ^ and 
b n = 6(°)fc~ 1/3 e- fc "+^ b ; here u(°) = 0.5, b^ = 0.05, and 
the random phases <p u and 4> b are distributed uniformly 
on the interval [— tt, ±7r]; different values of these random 
phases distinguish different initial conditions; we work 
with decaying turbulence, so there is no forcing term; 
and our results are averaged over 7500 independent 
initial conditions [TB]. We set N — 22, use a second- 
order, slaved Adams-Bashforth scheme [2U] for solving 
the shell-model ordinary differential equations ([2]), and 
calculate the energy spectra E u (k n ) = \\u n \ 2 /k n and 
E b {k n ) = \\b n \ 2 /k n (the superscripts u and b refer 
to velocity and magnetic field, respectively), the root- 
mean-square velocity u rms = \J^ n \ u n | 2 , the Taylor 
microscale A = ^jT, n E u {k n ) /T, n klE u (k n ) 7 the Taylor- 
microscale Reynolds number Re\ = u rms \/is c ff, the 
integral length scale Zj = £„ (E u (k n )/k n ) /H n E u {k n ), 
the effective viscosity and magnetic diffusivity 
y cS = E„ {vk 2 n E u (k n ) + y 2 k n E ll (k n )) /^ n k 2 n E u {k n ) 
and VcS = E„ (rjk 2 E b (k n ) + rj 2 k n E b {k n )) /E n fc 2 £ b (fc n ), 
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FIG. 1: (Color online) (a) Plots of e™ (red, upper curve) and e b (blue, lower curve) versus the rescaled time t/r. (b) Plot of the 
kinetic-energy spectrum E u (k); the thick, blue line indicates the K41-scaling power law. Inset : Time evolution of E u (k) during 
the initial period with transients; the legend indicates the time at which each curve is obtained, (c) Plot of the magnetic-energy 
spectrum E b (k); the thick, blue line indicates K41 scaling in the low-fc region and the dashed, black line indicates k~ i5 scaling 
in the high-fc region . Inset : Time evolution of E b (k) during the initial period with transients; the legend indicates the time 
at which each curve is obtained. Here we use E u (k n ) = ££(fc n )/fc„ and E\k n ) = T,\(k n )/k n to suppress 3-cycles (see text). 
The data shown are from run R2 and the spectral exponents for all our runs are reported in Table I similar plots for runs Rl, 
R3, and R4 are given in the Supplemental Material |21| . 



respectively, the effective magnetic Prandtl num- 
ber [55] PrMeS — ^cs/VcS, and the dissipation rates 
e u = v eS V n klE u (k n ) and e b = r, eS H n klE b (k n ) . The 
parameters of our simulations are given in Table I. 

In shell models, the equal-time, order-p structure func- 
tions for the velocity field and the magnetic field are 

defined, respectively, as S p (k n ) = (\u n \ p ) ~ k n " and 

S b p (k n ) = (\b n \P), where S b p (k n ) ~ k^ (k < kj) and 

Sp(k n ) ~ k n p (kd > k > ki). However, to remove the 
effects of an underlying three cycle in GOY-type shell 
models [TU [16] , we use the modified structure functions 
S p( fc n) = (\^bhi+2U n+1 u n + l/4u„_ 1 u„u„ +1 ]| p / 3 ) and 
E£(fc„) = (\%{b n+2 b n+1 b n + l/46 n _ 1 6 B & n+1 ]|i'/ 3 ), from 
which we can obtain multiscaling exponents via E^ (fc„) <~ 

kn , E^(fc n ) ^ k n p {kn ^ fciji and E^(fc n ) ^ k n p 
(kd > k n > kj). We also use the extended self-similarity 
(ESS) procedure [53] to determine exponent ratios from 
slopes of log-log plots of E^ versus Eg and their magnetic 
counterparts (Fig. 2 inset). 

In Fig. 1(a) we show plots of e u (red, upper curve) and 
e b (blue, lower curve) versus the rescaled time t/r, for run 
R2, where the box-size eddy-turnover time r = lj(u\k\) 
is evaluated at the principal peak of s u . This peak sig- 
nals the completion of the Richardson cascade 13], as we 
can see from the time evolution of E u (k n ) and E b (k n ), 
in the insets of Figs. 1(b) and (c), respectively, where the 
red lines with full circles denote the spectra at cascade 
completion. We evaluate the spectral-slope exponent a u 
at cascade completion from log- log plots of E u (k n ) ver- 
sus k n as shown in Fig. 1(b). We find a u ~ 5/3, as 
predicted by dimensional analysis, and illustrated in Fig. 
1(b) by a thick, blue line. In Fig. 1(c) we show a rep- 
resentative plot of E b (k n ); we see two different scaling 
regimes clearly: (1) from the low-fc one (solid, blue line), 
we find that a 6 ' 1 ~ 5/3, which is consistent with dimen- 
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FIG. 2: (Color online) The exponents £p (blue filled circles 
connected by a blue line), (^J' 1 (red filled squares connected by 
a red line), and the K41 prediction (thick black line) versus 
p. The lines connecting the data points are a guide to the 
eye. (Inset) A representative plot of EJJ versus £3; equal- 
time exponents are calculated from such plots. The different 
curves are for p — 1 (top) ... 6 (bottom) and are from run R2. 



sional analysis; (2) from the high-fc regime, we obtain 
a b ' 2 = 3.45 ± 0.06, which is close to the dimensional- 
analysis value 11/3 for IMHD systems, such as ours, in 
which induction and Hall terms are comparable |10L I12j . 
These spectral exponents are consistent with those in 
solar- wind experiments 7 j. In Table I we provide our 
results for all three spectral exponents; we obtain the 
values of these and all other exponents from the means 
of our runs with 7500 independent initial conditions; the 
error bars follow from the associated standard deviations. 

To characterize the statistical properties of the Hall- 
MHD system, we now calculate the equal-time exponents 
Cp, Cp 1 ( k < k i)i and Cp 2 {kd > k > ki) via the mod- 
ified structure functions S" and just after cascade 
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V 






a6,2 
S>p 


At., 2 i r b,1 
S>P / t >3 


1 


0.36 ±0.01 


0.348 ±0.006 


1.40 ±0.04 


0.38 ±0.01 


2 


0.69 ±0.01 


0.682 ±0.006 


2.6 ±0.1 


0.698 ±0.009 


3 


1.0 


1.0 


3.7 ±0.2 


1.0 


4 


1.29 ±0.02 


1.31 ± 0.01 


4.8 ±0.3 


1.30 ±0.01 


5 


1.57 ±0.06 


1.60 ±0.03 


6.0 ±0.4 


1.60 ±0.03 


6 


1.85 ±0.09 


1.89 ±0.05 


7.1 ±0.4 


1.90 ±0.04 


7 


2.1 ± 0.1 


2.18 ±0.08 


8.2 ±0.6 


2.20 ±0.06 


8 


2.4 ±0.2 


2.5 ±0.1 


9.3 ±0.6 


2.50 ±0.07 


9 


2.7 ±0.2 


2.7 ± 0.1 


10.5 ±0.7 


2.81 ±0.08 


10 


3.0 ±0.3 


3.0 ±0.2 


11.6±0.8 


3.11 ±0.09 



TABLE II: Multiscaling exponents Q, Q' 1 , Q> 2 , and the ex- 
ponent ratio Cp' 2 /C3' 2 f° r run R-2. The exponents from the 
other runs are equal to the ones shown here. We note that 

q ~ Cp 1 ' 1 * C^'Vcl'V & 2 ED- 




FIG. 3: (Color online) Plot of the exponents Cp' 2 versus the 
order p showing a mild departure from linear scaling. Inset 
: the exponent ratios Cp'VCs' 2 (obtained via ESS) versus p, 
showing clear multiscaling. The thick, black line indicates 
simple scaling. The Exponents Cp' 2 and Cp'VCs' 2 are shown 
for runs R1-R4 (see legend). 



completion. We find C3 1 = 1 = Q 1 , which is consis- 
tent with dimensional analysis. In Fig. 2 we show the 
order-p equal-time exponents Cp (blue, filled circles) and 
Cp' 1 (red, filled squares) for integer values of p between 1 
and 10; the thick, black line illustrates the dimensional, 
simple K41 scaling. We see that Cp ~ Cp' 1 and both expo- 
nents show clear multiscaling corrections to K41 scaling, 
with values consistent with those obtained in 3D MHD 
turbulence [3]. We obtain these multiscaling exponents 
by using Ep 1 and Ep and the ESS procedure [23], to ex- 
tend the scaling range. However, the result C3 ~ 1 en- 
sures that the exponent ratios and the exponents them- 
selves are equal (within error bars). In Table II, we list 
the order-p equal-time exponents Cp and C, for integer 
values of p between 1 and 10. 

We finally turn to the exponents Cp' 2 , which character- 
ize the high- k regime (fc<j > k > kj). Solar-wind measure- 
ments [7] of Cp' 2 suggest simple-scaling behaviour, with 
Cp' 2 a linear function of p. In Fig. 3 we plot Cp 2 (ob- 
tained without ESS) versus p. Our results for Cp 2 are 
in qualitative agreement with solar-wind measurements 
to the extent that there is only mild multiscaling; i.e., 
Cp' 2 is a nonlinear, monotone, increasing function of p, 
but the deviation from a linear dependence on p is not 
very pronounced. Although the exponents Cp' 2 , for all 
the runs R1-R4, are in agreement with each other, given 
our error-bars, their mean values seem to decrease with 
dj. We now use the ESS procedure to obtain the expo- 
nent ratios Cp'VCs' 2 , which are plotted versus p in the 
inset of Fig. 3 (note C3' 2 7^ 1). This ESS plot is remark- 
able for two reasons : (1) There is a clear signature of 
multiscaling (the thick, black line in the inset indicates 
simple scaling); (2) although the exponents Cp' 2 are very 
different from Cp and Cp \ the ratios Cp VC3' 2 are equal 
to Cp and Cp' 1 (within error bars). In Table II, we list 



both Cp 2 and CpVCss' 2 f° r different values of p for the 
representative run R2. 

We hope our extensive studies of the multiscaling of 
structure functions in a shell model for 3D Hall-MHD will 
stimulate high-precision and high-resolution experimen- 
tal and DNS studies to determine conclusively whether 
3D Hall-MHD turbulence shows multiscaling for kd > 
k > kj. Our ESS results suggest that structure func- 
tions show mild, but distinct, multiscaling in this region 
To obtain quantitative agreement with solar-wind expo- 
nents, we must, of course, carry out DNS studies of the 
3D Hall-MHD equations ([T]) and include compressibility 
effects and a mean magnetic field [SHlOj; however, current 
computational resources limit severely the spatial resolu- 
tion of such DNS studies so they cannot (a) uncover the 
multiscaling of magnetic-field structure functions in 3D 
Hall-MHD turbulence in both low- and high-/c power-law 
ranges and (b) obtain well-averaged multiscaling expo- 
nent ratios. For the moment, therefore, the shell-model 
study, which we have undertaken, provides the only way 
of understanding the multiscaling of structure functions 
in the solar wind [7] and the apparent and intriguing 
universality of the exponent ratios. This apparant uni- 
versality needs to be investigated in detail in experiments 
and DNS. 
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SUPPLEMENTAL MATERIAL 

In this Supplemental Material we describe details of 
our work that are of interest only to specialists of the 
field. The notations and abbreviations used in this Sup- 
plemental Material are the same as in the main paper. 

In our main paper, we discuss the statistical nature of 
turbulence in the Hall-MHD plasma and, in particular, 
the scaling properties of various structure functions, in 
great detail. To substantiate our claims, we show in the 
main paper representative data from only a single set of 
simulations (except in Fig. 3 where we show results from 
all our simulations), namely, Run R2 (see Table I of the 
main paper). In Fig. 3 of the main paper we do show 
exponents from all the four different sets of simulation 
(as detailed in Table I of the main paper); however, in 
Table II (main paper) we list exponents from R2 only; 
because the different sets of simulations all agree with 
each other, representative data from one set of simula- 
tions (R2), in the main paper, is enough to highlight the 
nature of multiscaling in Hall-MHD turbulence. 

We give here the equal-time exponents from the runs 
Rl, R3, R4 in Table I; these exponents are in agreement 
with the ones listed for R2 in Table II of the main paper. 
Furthermore, we show plots for the fluid and magnetic 
energy dissipation rates (Fig 1(a) for run Rl, Fig 2(a) 
for run R3, and Fig 3(a) for run R4), the kinetic energy 
spectrum (Fig 1(b) for run Rl, Fig 2(b) for run R3, and 
Fig 3(b) for run R4), and the magnetic energy spectrum 
(Fig 1(c) for run Rl, Fig 2(c) for run R3, and Fig 3(c) 
for run R4); these plots are analogous to the plots shown 
in Fig.l, for run R2, in the main paper. 

The parameters of the runs Rl, R3, and R4 (along 
with those for run R2) are given in Table (I) of our main 
paper. 



G 





Rl 


R3 


R4 


p 


Q/G 




/-b,2 i r b,2 
V / ^3 


q/q 




^P 1 ^3 






S>p /S3 


1 


0.36 ±0.01 


0.349 ± 0.007 


0.38 ±0.01 


0.36 ±0.01 


0.348 ±0.005 


0.37 ±0.01 


0.355 ±0.007 


0.349 ± 0.005 


0.38 ±0.01 


2 


0.69 ±0.01 


0.682 ± .007 


0.70 ± .01 


0.69 ±0.01 


0.681 ± .005 


0.69 ± .01 


0.687 ±0.008 


0.682 ± .005 


0.70 ± .01 


3 


1 


1 


1 


1 


1 


1 


1 


1 


1 


4 


1.29 ± 0.02 


1.30 ±0.01 


1.30 ±0.01 


1.30 ±0.02 


1.306 ±0.009 


1.30 ±0.01 


1.30 ±0.02 


1.31 ±0.01 


1.30 ±0.01 


5 


1.57 ±0.06 


1.60 ±0.03 


1.61 ±0.03 


1.58 ±0.04 


1.60 ±0.02 


1.61 ±.03 


1.58 ±0.04 


1.60 ±0.03 


1.60 ± .03 


6 


1.8 ±0.1 


1.89 ±0.05 


1.91 ±0.04 


1.87 ±0.07 


1.89 ±0.04 


1.91 ±0.05 


1.87 ±0.06 


1.89 ±0.04 


1.91 ±0.05 


7 


2.1 ± 0.1 


2.18 ±0.08 


2.22 ±0.06 


2.1 ±0.1 


2.18 ±0.05 


2.22 ±0.07 


2.15 ±0.08 


2.18 ±0.07 


2.22 ±0.06 


8 


2.4 ± 0.2 


2.4 ±0.1 


2.53 ±0.07 


2.4 ±0.1 


2.46 ±0.07 


2.52 ±0.09 


2.4 ±0.1 


2.47 ±0.09 


2.52 ±0.08 


9 


2.7 ± 0.2 


2.7 ±0.1 


2.84 ±0.09 


2.7 ± 0.2 


2.74 ± 0.09 


2.8 ±0.1 


2.7 ± 0.1 


2.8 ±0.1 


2.8 ±0.1 


10 


3.0 ±0.2 


3.0 ±0.1 


3.2 ±0.1 


3.0 ±0.2 


3.0±0.1 


3.1 ±0.1 


3.0 ±0.2 


3.0 ±0.1 


3.1 ±0.1 



TABLE III: Multiscaling exponent ratios Cp/Cs, Cp'Vd' 1 , and Cp 2 /Cs' 2 with error bars for the order p in the range 1 < p < 10 
and runs Rl, R3, and R4. 




FIG. 4: (Color online) (a) Plots of e u (red, upper curve) and e b (blue, lower curve) versus the rescaled time t/r. (b) Plot of the 
kinetic-energy spectrum E u (k); the thick, blue line indicates the K41-scaling power law. Inset : Time evolution of E u (k) during 
the initial period with transients; the legend indicates the time at which each curve is obtained, (c) Plot of the magnetic-energy 
spectrum E b (k); the thick, blue line indicates K41 scaling in the low-fc region and the dashed, black line indicates fc -3,5 scaling 
in the high-fc region . Inset : Time evolution of E b (k) during the initial period with transients; the legend indicates the time 
at which each curve is obtained. The data shown are for run Rl. 




FIG. 5: (Color online) (a) Plots of e™ (red, upper curve) and e b (blue, lower curve) versus the rescaled time t/r. (b) Plot of the 
kinetic-energy spectrum E u (k); the thick, blue line indicates the K41-scaling power law. Inset : Time evolution of E u (k) during 
the initial period with transients; the legend indicates the time at which each curve is obtained, (c) Plot of the magnetic-energy 
spectrum E b (k); the thick, blue line indicates K41 scaling in the low-fc region and the dashed, black line indicates fc -3 ' 5 scaling 
in the high-fc region . Inset : Time evolution of E b (k) during the initial period with transients; the legend indicates the time 
at which each curve is obtained. The data shown are for run R3. 
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FIG. 6: (Color online) (a) Plots of e u (red, upper curve) and e b (blue, lower curve) versus the rescaled time t/r. (b) Plot of the 
kinetic-energy spectrum E u (k); the thick, blue line indicates the K41-scaling power law. Inset : Time evolution of E u (k) during 
the initial period with transients; the legend indicates the time at which each curve is obtained, (c) Plot of the magnetic-energy 
spectrum E b (k); the thick, blue line indicates K41 scaling in the low-fc region and the dashed, black line indicates k~ i5 scaling 
in the high-fc region . Inset : Time evolution of E b (k) during the initial period with transients; the legend indicates the time 
at which each curve is obtained. The data shown are for run R4. 



